Time-Resolved Magnetic Relaxation of a Nanomagnet on Subnanosecond Time Scales 
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We present a two-current-pulse temporal correlation experiment to study the intrinsic subnanosec- 
ond nonequilibrium magnetic dynamics of a nanomagnet during and following a pulse excitation. 
This method is applied to a model spin-transfer system, a spin valve nanopillar with perpendicular 
magnetic anisotropy. Two-pulses separated by a short delay (< 500 ps) are shown to lead to the 
same switching probability as a single pulse with a duration that depends on the delay. This demon- 
strates a remarkable symmetry between magnetic excitation and relaxation and provides a direct 
measurement of the magnetic relaxation time. The results are consistent with a simple finite tem- 
perature Fokker-Planck macrospin model of the dynamics, suggesting more coherent magnetization 
dynamics in this short time nonequilibrium limit than near equilibrium. 



The control of magnetization on short timescales has 
become an area of intense research and many methods 
are used to excite a magnetic system, including spin- 
currents [l|, [H , magnetic fields |3|] and optical pulses 0, H| . 
It also has recently become possible to drive individual 
nanometer scale magnetic elements far from equilibrium 
using spin-currents and probe their dynamical response. 
However, most electronic transport studies focus on the 
excitation during the drive pulses @, @] • Less is known 
about the relaxation processes following a pulse excita- 
tion. Yet, this is very important for fundamental and 
practical reasons. First, the time scales of the relaxation 
far from equilibrium are not known and may differ from 
those determined in experiments that probe low ampli- 
tude magnetic excitations, such as ferromagnetic reso- 
nance. Second, these time scales determine the speed 
at which nanomagnets can be written and read in mag- 
netic memories, including in spin- transfer torque mag- 
netic random access memories (STT-MRAM). 

It is now well known that spin-polarized currents and 
spin-transfer torques (STTs) can reverse the direction of 
the magnetization on subnanosecond timescales [H, Q. 
In fact, many experimental methods have been devel- 
oped to study spin-transfer switching in both spin valves 
and magnetic tunnel junctions (MTJs). While MTJs are 
promising for applications because of their large magne- 
toresistive (MR) readout signals (> 100%), all- metallic 
spin valves permit one to apply much larger currents 
and thus can be driven farther away from equilibrium. 
Direct single-shot time-resolved electrical measurements 
have been carried out in MTJs @, 0] , but have not been 
reported in spin valve nanopillar s, because of their low 
impedance and small magnetoresistance (MR) (< 5%), 
which results in their switching signals being to small for 
nanosecond measurements. Further, time-resolved STT 
measurements are insensitive to the magnetization dy- 
namics after the excitation, as the excitation current is 



the source of the electrical readout signal. Therefore, to 
time-resolve the relaxation dynamics, a different experi- 
mental method is needed. 

In this paper we present an all-electrical two-pulse cor- 
relation method that yields quantitative information on 
the form and timescales of magnetic excitation and re- 
laxation, with 50 ps time resolution. This method is 
applied to a model spin-transfer system, a spin valve 
nanopillar with perpendicular magnetic anisotropy. Two- 
pulses separated by a short delay (< 500 ps) are shown 
to lead to the same switching probability as a sin- 
gle pulse with a duration that depends on the delay- 
demonstrating a remarkable symmetry between magnetic 
excitation and relaxation-and providing a direct mea- 
surement of the magnetic relaxation time. The observed 
symmetry is shown to be consistent with a finite temper- 
ature macrospin model. 

The system we studied has two stable magnetic states 
A and £?, which are nearly degenerate in equilibrium. We 
are interested in how the system relaxes to these states 
after being excited by a spin-polarized current. The basic 
idea of our method is to compare the switching proba- 
bility of two pulses separated by a delay much longer 
than the magnetic relaxation time to that of a composite 
pulse, i.e. two pulses separated by a short delay. With a 
composite pulse, the first pulse alters the magnetic state 
on a time scale that changes the dynamics excited by the 
second pulse. By choosing the amplitudes and polarities 
of the first and second pulse we can study different as- 
pects of the relaxation processes. For example, to study 
the relaxation back to the initial state (i.e. state A) af- 
ter a pulse, we can apply two pulses with identical po- 
larities that both would be of the appropriate polarity 
to switch the sample from state A to B. If the delay 
is much longer than the relaxation time £ re iax, the two 
pulses can be considered to be independent events and 
the total probability of switching from state A to state 
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FIG. 1: (a) The total switching probability, Pap^+p, as a function of the second pulse duration £2, for a current I\ of 8.5 mA 
and different delays. The dashed line is the switching probability P^ P ^ P for a single 8.5 mA pulse of duration £2. (b) The 
same data is plotted with shifted time axis, with the shift dependent on the delay ££(£deiay)- The data closely overlaps showing 
that the switching probability of a double pulse (£1 + £dei ay + £2) is the same as that of a single pulse with a longer duration, 
i.e. £2 + ££(£deia y ), as illustrated schematically by the pulses drawn within the figure. 
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Here P^\ B ^ P^Xb anc ^ Pa^*b g are the switching proba- 
bilities from state A to B for the first, the second and for 
both pulses combined and we have assumed no reverse 
switching during the measurement, i.e. Pb^a = 0. If 
^deiay is within the relaxation time of the system, the total 
switching probability will increase compared to Eq. (j2j): 
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Using these equations we can resolve the relaxation dy- 
namics of the sample by measuring and comparing the 
switching probabilities as a function of £deiay [S|- 

Our sample is a spin-valve nanopillar consisting of a 
Ni|Co free layer and exchange coupled Ni|Co and Co|Pt 
multilayers as the reference layer, both of which have 
perpendicular magnetization anisotropy (PMA). (For the 
full layer stack see 0, Hoj|.) Such PMA samples have a 
simple uniaxial anisotropy energy landscape that permits 
a direct comparison to theoretical models. The sample 
presented in this work is 100 nm x 100 nm in size, fab- 
ricated by e-beam and optical lithography 

We 

have studied more that 20 samples of different shapes and 
sizes and found similar results. The equilibrium states 
are antiparallel (A = AP) and parallel (B = P) magneti- 
zation alignment of the free and reference layers and there 
is 0.3% MR, well within the accuracy of our method. The 
room temperature coercive field is 100 mT, large enough 
so that no reverse switching occurs [13] . The zero temper- 
ature switching current for antiparallel (AP) to parallel 
(P) switching is I c q = 6.5 mA. 



We use an arbitrary waveform generator to apply two 
pulses separated by a time delay tdeiay, as illustrated by 
the inset of Fig. QJa). The state of the spin valve is 
determined by a resistance measurement using a lock- 
in amplifier both before and after the injection of the 
double-current pulse. We note that the lock-in current is 
small enough (300 /iA) not to affect the spin valves and 
not to induce reverse switching jTjj. All experiments 
described here were conducted at room temperature. 

Initially the sample was brought into the AP state by 
the application of an easy-axis magnetic field of appro- 
priate direction and magnitude. A short current pulse of 
positive polarity I\ = 8.5 mA and t\ = 0.9 ns was then 
applied. A longer pulse of positive polarity would even- 
tually switch the sample into the P state (as I\ > I c q). 
However, the duration has been chosen such that the 
switching probability is nearly zero (< 1/100), as shown 
by the dashed curve (i.e. the single pulse curve) in Fig. 
QJa) at 0.9 ns. After a variable delay from 50 ps to 1.5 ns, 
a second pulse with the same amplitude (I2 = h) and a 
variable duration £2, from 50 ps to 5 ns, was applied. The 
total switching probability distribution P^p^p was then 
measured by repeating the process 100 to 10 000 times 
for each pulse set. 

As seen in Fig. [U^a), for delays longer than 500 ps the 
switching probability distribution is the same as that of 
the second pulse alone P^^ e = P^p^p: which is con- 
sistent with Eq. (j2]) for P^\b = 0- This indicates that 
these delays are longer than the magnetization relaxation 
time. The magnetization state excited by a 0.9 ns cur- 
rent pulse decays to equilibrium in under 500 ps, within 
our experimental accuracy. However, for delays less than 
500 ps, the switching probability P^p^p is larger than 
that of just a single pulse. 
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Remarkably, even with a composite pulse, the switch- 
ing probability versus the second pulse duration is seen 
to follow the same functional form as that for a single 
pulse. This is shown in Fig. QJb), where the probabil- 
ity data, shifted by an amount St (t delay) ~ a quantity 
that depends on the delay between the pulses - is seen 
to closely overlap. This clearly demonstrates that the 
switching probability of a composite pulse (t\ + tdeiay +^2) 
is equal to that of a single pulse with an increased dura- 
tion (t2 + St). Since this is satisfied for all £2, a natural 
conclusion is that the ensemble averaged state of the sam- 
ple after the first pulse and the delay (t\ + tdeiay) is the 
same as that after a single pulse with a duration of St, as 
illustrated schematically by the pulse shapes in the inset 
of Fig.lHb). 

From the data shown in Fig. [1] we determine St as a 
function of the delay f delay This is shown as crosses in 
Fig. El The black solid curve is an exponential fit to the 
data where tl is the lifetime of the excitation: 



St = to exp 



^delay 
TL 



(3) 



with the parameters to = 0.95 ns and tl = 0.28 ns. For 
no delay, tdeiay — 0, clearly St = t\ and thus to = ti, 
which is indeed fulfilled within the time resolution of the 
experiment (p± 50 ps). Eq. (j3j) links the duration of the 
excitation, the delay and the corresponding equivalent 
pulse duration, which is related to the magnetization re- 
laxation rate. 
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FIG. 2: St versus delay tdeiay The blue crosses are deter- 
mined from the experimental data in Fig. [1] The black solid 
curve is an exponential fit to Eq. ©. The green dash-dotted 
line is a calculation from the LLG equation for T = 0. The 
red dashed curve is a calculation based on a Fokker-Planck 
Analysis (Eq. ©). 

The data in Fig. [2] is unexpected at first, as a zero 
temperature macrospin Landau-Lifshitz-Gilbert (LLG) 
model predicts a linear relationship between St and tdeiay, 
as shown by the green dash-dotted line in Fig. O This 
linear relationship follows from the LLG equation. For a 
uniaxial nanomagnet in zero applied magnetic field, the 



polar angle of the magnetization i.e. the angle between 
the magnetization and the easy axis, satisfies [lij : 



dO 
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T£> is the relaxation time due to damping and depends 
on the material parameters, the anisotropy field H^, the 
Gilbert damping parameter a and the gyromagnetic ra- 
tio 7. td is typically hundreds of picoseconds. I c o is the 
zero temperature critical current as we mentioned ear- 
lier, representing the threshold current for switching. It 
depends on the magnetization density M s , the nanomag- 
net's volume V and the spin polarization of the current 
77. From Eq. both excitation (I > I c o) and relaxation 
(I = 0) follow the same functional form for small angles 
and are only different in their respective time constants. 
It is straightforward to show that the duration of a single 
pulse (St), which brings the sample to the same state as 
the first pulse (to) followed by the delay (tdeiay) satisfies: 



St = to 



IcO 



^delay 
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Hence the initial slope of St versus tdeiay? shown as the 
green dashed dotted curve in Fig. [2j depends on the cur- 
rent overdrive, I/I c o ~ 1? an d provides an independent 
method to determine the zero temperature critical cur- 
rent. However, for long delays this equation is unphysi- 
cal, as it predicts a negative St. 

A physical picture is that during the first pulse the 
magnetization evolves from a thermal distribution of ini- 
tial angles near the north pole ((0 O ) ~ y/irk B T/W « 6°, 
where U is the energy barrier and /c#T is the thermal 
energy) to a larger angle 6\. During the delay, and pro- 
vided 0i < 90°, the polar angle decays back toward the 
north pole. However, in this zero temperature model the 
magnetization decays back to O2 = 0, i.e. to an angle 
that can be less than the initial angle, #o- 

The difference between the LLG model prediction and 
the experimental results can be explained by thermal 
fluctuations both during the excitation and relaxation 
processes. As is well known, thermal fluctuations lead 
to a spread of the initial angles of magnetization about 
the north pole. A spin-current pulse works to rotate the 
magnetization away from the easy axis direction and, in 
essence, amplifies these magnetization fluctuations. Once 
the magnetization deviates from the easy axis direction, 
the magnetization polar angle grows exponentially. How- 
ever, thermal noise also slows down the relaxation from a 




FIG. 3: Results of a Fokker-Planck calculation, (a) Ensemble averaged magnetization projection versus time for a 1 ns delay. 
During the first pulse the average magnetization projection increases with time as shown by the first red curve, and it decreases 
with time during the delay as shown by the blue curve. The shaded region depicts the angular range containing 10% to 90% 
of the states, (b) Occupation probability plotted as a function of the polar angle. The green curve is the initial thermal 
equilibrium Boltzmann distributed polar angle. The red solid curve is the probability distribution at the end of the first pulse. 
The light blue solid curve is the probability distribution after the first pulse and delay, tdeiay = 0.1 and 0.5 ns. The black 
dashed curve is the probability distribution corresponding to a pulse of duration St. The corresponding solid light blue and 
black dashed curves match well showing that these distributions are nearly same - that the probability distribution after a 
pulse and delay is nearly the same as that after a single shorter pulse St - as inferred from the scaling of the experimental data 
presented in Fig. 1(b). (c) Switching probability distributions plotted as a function of the second pulse duration for different 
delay times. As in the experiments, the probability distributions have the same form but are shifted horizontally, along the 
time axis. 



simple exponential decay process predicted by the LLG 
equations and leads to relaxation to a thermal distribu- 
tion. Therefore, thermal fluctuations not only change 
the timescales of the excitation and relaxation processes, 
they also alter their functional forms. 

To model the influence of thermal fluctuations, we sim- 
ulated the double pulse measurement using the Fokker- 
Planck Equation [15j , which has the following form for a 
uniaxial macrospin system: 



dV 
~dt 
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where V = V (0, t) is the probability density of the mag- 
netization at the angle and the time t. £ is the di- 
mensionless energy barrier, £ = [//(fc^T). For a uniaxial 
nanomagnet U = l/2M s HkcV . 

We start with a Boltzmann distribution centered 
around the easy axis in the AP state: 



p(e,t = o) 



Po exp (£ cos 2 6) sin 




I 



0° < 6 < 90° 
90° < < 180° 



Where Po is a normalization constant so that 

V (0, t = 0) dO = 1. For the simulation we use a fixed 

duration for the first pulse of 0.9 ns, the same as in our 
experiments and vary the duration of the second pulse 
for different delays. By fitting the measured single pulse 
switching probability distribution (i.e. the red dashed 
curve of Fig. E{a)) to Eq. (|8j). we obtain £ = 450 and 
T£> = 354 ps. 



We then can track the ensemble averaged magnetiza- 

r 

tion projection, (cosO) = / cosOV (6,t) d6 as a func- 

Jo 

tion of time f, as shown in Fig. [3fa). The shaded region 
depicts the angular range containing 10% to 90% of the 
states. During the first current pulse (0 ns < t < 0.9 ns), 
the switching process begins. The averaged magnetiza- 
tion projection decreases while the occupation probabil- 
ity distribution spreads out, where the majority of the 
distribution is still centered around = 0°. During the 
delay (0.9 ns < t < 1.9 ns), the averaged magnetization 
projection increases while the occupation probability dis- 
tribution narrows toward the initial state. Then during 
the second pulse (t > 1.9 ns), the averaged magnetiza- 
tion projection decreases again as the total distribution 
expands during the beginning part of the switching pro- 
cess. 

We can further check whether a pulse and a delay 
(ti + ^deiay) brings the system into the same state as 
a single St duration pulse. Fig. [3fb) shows the occu- 
pation probability as a function of the polar angle for 
different delay times. The distributions for a pulse and 
delay t\ + tdeiay (the solid light blue curves in Fig. EJb)) 
overlap with those of a single St pulse (dashed curves in 
Fig. [3jb))-which is consistent with scaling of the experi- 
mental data shown in Fig. QJb). 

The Fokker-Planck simulation results also agree quan- 
titatively with the experimental data. As shown in 
Fig. [3fc) the calculated switching probability distribu- 
tions have the same form independent of tdeiay and can 
be shifted in time to coincide. The results shifts, Sts, 
are in good agreement with the experimental results, as 
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indicated by the red dashed curve in Fig. [2j The differ- 
ence between the LLG and the Fokker-Planck analysis 
clearly shows the influence of the thermal fluctuations, 
even though the total process (excitation plus relaxation) 
lasts less than 2 ns. 

In summary, we have presented a method ideally suited 
to study excitation and relaxation in a magnetic nanos- 
tructure. This method only requires the generation of 
pulses with short durations and delays and measurement 
of the switching probability. It does not require resolving 
small signal changes at high bandwidth and it can also 
be used to study dynamics after a pulse, with no current 
applied. 

Our method provides a direct (model independent) 
measurement of the relaxation timescale by consider- 
ing when two events are correlated. We have used this 
method to study excitation and relaxation in a model 
system-a thin film nanomagnet with uniaxial anisotropy 
excited by spin-current pulses. We have shown that a 
single pulse followed by a delay leaves a nanomagnet in 
a state with the same ensemble average as that of a sin- 
gle pulse with a shorter duration St. This finding is also 
model-independent. The Sts are seen to decrease expo- 
nentially, following the form of Eq. [3j and giving a relax- 
ation time of 280 ps. We note that the relaxation time is 
more typically obtained experimentally by fitting to an 
Arrhenius-Neel's model with a variable prefactor (see, for 
example, [HI). ^ n general, the decay of the magnetization 
might be expected to be more complex than a simple ex- 
ponential (Eq.[3)), including multiple relaxation pathways 
with different timescales. However, even in such a case 
the directly measured relationship between St and tdelay 
would provide information that is important for under- 
standing such relaxation processes. Using Eq. (j5j) with 
an anisotropy field of Hk = 0.25 T (see Ref. ^) and a 
gyromagnetic ratio of 7 = 1.76 x 10 11 /(Ts) (i.e. a Lande 
g- facto of 2), we find a Gilbert damping of a = 0.09, a 
value that is larger than that determined from thin film 
ferromagnetic resonance measurements of the damping 
(a = 0.04 (l7j). The larger damping inferred from the 



data may be a consequence of the large-angle excitation 
of the magnetization, including non-linear effects [HI, or 
to changes of the magnet's material properties associated 
with device nanofabrication [l9| . 

Interestingly, our results can be understood within a 
simple macrospin model that incorporates thermal fluc- 
tuations. In fitting our data to this model we obtain a 
dimensionless energy barrier of £ = 450, which is within 
25% of that expected based on the volume, magnetization 
and anisotropy of the nanomagnet (£ = M s HkV/ (2fcT) = 
360). Our previous results, as well as those of many 
other groups, have found energy barriers to reversal for 
long field [2CH23j and long current pulses (> 100 ns) 
[l, [loj to be only a small fraction of that expected for 
a macrospin, sug gest in g the reversal proceeds by subvol- 
ume nucleation 24], 25j. For similar samples, we found 



that the energy barrier to reversal under long current 
pulses to be ~ 60 keT [9], 1/6 the macrospin value. 
This suggests that the reversal modes in the short-time 
nonequilibrium limit are distinct from those in the long- 
time limit, as the system may not have time to follow 
minimum energy paths or be driven from these paths by 
the spin transfer torque. The fundamental explanation 
is an open question. We expect that the dependency 
of the effective energy barrier on pulse amplitude can be 
further studied by our method with different pulse condi- 
tions. Our double pulse method, which provides access to 
the intrinsic time scales of the excitation and relaxation 
dynamics, can be used in other cases where the readout 
signal is also small, such as current induced domain wall 
motion. 
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